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Abstract. We prove a version of differential Harnack inequality 
for a family of sub-elliptic diffusions on Sasakian manifolds under 
certain curvature conditions. 

1. Introduction 

Harnack inequality is one of the most fundamental results in the the- 
ory of elliptic and parabolic equations. For linear parabolic equations 
in divergence form, this was first done in p3j. Since then, numerous 
developments around this inequality were found. In [13], the, so called, 
Li-Yau estimate was proved. This is a sharp gradient estimate for linear 
parabolic equations on Riemannian manifolds with a lower bound on 
the Ricci curvature. This estimate is also called a differential Harnack 
inequality since one can recover the Harnack inequality by integrating 
this estimate along geodesies. 

There are many generalizations of the Li-Yau estimate for geometric 
evolution equations. This includes the evolution equations for hyper- 
surfaces [HJ [7J [T], the Yamabe flow [8], the Ricci flow [10] and its 
Kahler analogue [5]. For a more detail account of these generalizations 
as well as further developments, see [15] . 

There are also generalizations [2] of the Li-Yau estimate to linear 
parabolic equations of the form 

(1.1) Pt = LfH 

under certain conditions called curvature-dimension conditions. Here 
L is a second order linear elliptic operator without constant term. The 
curvature-dimension conditions were recently generalized by [3] to ob- 
tain Li-Yau type estimates for equations of the form (11.11) . where L is 
a linear sub-elliptic operator without constant term. The following is 
one of the main results in [3] when the underlying manifold is Sasakian 
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and the equation is the sub-elliptic heat equation (for the definition of 
Sasakian manifolds and various related notions appeared in Theorem 
HH see Section 2). 



Theorem 1.1. [3] Assume that the manifold is Sasakian and satisfies 
Rc > 0. Then any positive solution of the sub- elliptic heat equation 

p t = A hor p t 

satisfies 

i + -}f t + l\v h0 M 2 +^\v ver f t (x)\ 2 < 2n( \ + - )2 



nj 2' 1 3 

for all t > and all x in M, where f t = — 21ogp 4 . 

We remark that a local version of the above estimate, which did not 
take the curvature into account, appeared in an earlier work [5|. 

On the other hand, the author proved in [12] another version of the 
differential Harnack inequality which is different from that of [2J. 

Theorem 1.2. [12] Let pt be a positive solution of the equation 

Pt = A Pt + (Vp t , Wi) + U 2 p t 

on a compact Riemannian manifold of non-negative Ricci curvature. 
Assume that 

A f-AUi - ^| Wi| 2 + 2U 2 ) > k 3 . 



Then 
where 




2A\ogp t + AU 1 > -na M (t), 

71 

ifK>0 

<f ifK = 

V^K coth( v /z X t) if K < 0. 

In this paper, we combine the ideas from [3] and [12] to obtain a 
differential Harnack estimate for fll. ID . where L is a linear sub-elliptic 
operator (possibly with constant term). In fact, we allow L to have a 
mild non-linearity (see Theorem 12.11 for the detail). In the case when 
the manifold is Sasakian and L is linear, we have the following result. 

Theorem 1.3. Assume that the manifold is a compact Sasakian and 
it satisfies Rc > 0. Let U\ and U2 be two smooth functions on M 
satisfying 

(1) V = A hor U 1 + \\V hor U l \ 2 -2U 2 

(2) V<K lf 
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(3) A hor V + £ (1 + f ) 2 1 V ver V\ 2 < k 2 , 

for some positive constants K\ and n 2 . Then any positive solution of 
the equation 

pt = A hor p t + (VhorUl, VhorPt) + U 2 p t 

satisfies 

^ ft{x) + \\Vkorft\l-V{x) 



n J ' T 



+ f 1 + Ti if^t a nh(c 2 t)\\/ ver f t (x)\ 2 < — coth(c 2 t) + — 
V 3/ y 2k 2 c 2 n 

for allt >0 and all x in M, where c 2 = ^3 V^lF an< ^ f f = — ^ A — 

As usual, we can integrate the estimate in Theorem 11.31 and obtain 
a Harnack estimate. Let 

II7 , . / n \ 2 Jr , s 3nK 4 

W(x) = V(x) H V 

K J \n + 3j K J (n + 3) 2 

In this case, we consider the following cost function. 

rsi y 

c S0>Sl (x ,xi) =inf / -|7(s)| 2 + W(~i(s))ds 

where / ranges over all smooth curves 7(-) such that j(s) G -D 7 ( s ) for 
all t in [sq, si]. 

Corollary 1.4. Assume that the manifold is compact Sasakian and 
it satisfies Rc > 0. Let JJ\ and U 2 be two smooth functions on M 
satisfying 

(1) V = A hor U 1 + ±\V hor U 1 \ 2 -2U 2 

(2) V < ki, 

(3) A hor V + £ (1 + 1 ) 2 |v„ er y| 2 < k 2 , 

for some positive constants K\ and n 2 . Then any positive solution of 
the equation 

pt = A hor p t + (VhorUl, VhorPt) + U 2 p t 

satisfies 



Ps 1 (x 1 ) ( sinh(c 2 s 



j \ ~(n+3) 

Ps (xo) ~ \smh(c 2 s ) J 

exp ( - i ^(xO - tf(z„) + (l + ^ 
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for all s , Si > and all x , Xi in M, where c 2 = ^3 y^F- 

The structure of the paper is as follows. The main results of this 
paper are stated in Section 2. In [12], a moving frame argument was 
used for the proof of Theorem 11.21 instead of the Bochner formula. 
The advantage is that a matrix version of Theorem ll.2[ generalizing 
the matrix Hamilton-Li- Yau estimate for the heat equation [9], can be 
proved using a very similar argument. Although there is no matrix 
analogue of Theorem 11.31 in this paper, we show that a version of the 
moving frame argument is possible in the present setting. This is done 
in Section 3 and 4. The proofs of the main results are given in Section 
5. Section 6 is an appendix devoted to some calculations needed in the 
proofs. 

2. The main results 

In this section, we give the statements of the main results. First, let 
us introduce the setup which is essentially the same as that of [3]. 

Let M be a Riemannian manifold and let us fix a distribution D (a 
vector bundle of the tangent bundle TM) of rank k. We assume that 
the orthogonal complement of D is spanned by n — k vector fields de- 
noted by wi, w n _fc which satisfy certain symmetry conditions to be 
specified. We will also assume that the vector field X t is the horizontal 
gradient Vh or /i of a one-parameter family of functions f t defined on 
the manifold M and specialize Lemma [3.41 to this case. 

Let us call vectors or vector fields which are contained in the distri- 
bution D horizontal. Let tpt be the flow of a vector field w. Assume 
that ipt sends horizontal vector fields to horizontal ones and preserves 
their lengths. If X\ and X 2 are horizontal vector fields, then we have 

((Vt)*(Xx),v) =0 and (^)*(X 2 )> = (Xi,X 2 ). 

for any vector field v which is in the orthogonal complement of D. 

If we differentiate the above equations with respect to t, then we 
obtain 

(\w,X 1 ],w) = and (V Xl w,X 2 ) + (X u V Xa w) = 0. 

Therefore, we call a vector field w which satisfies the following two 
conditions horizontal isometry: 

• [wj,Xi] is horizontal, 

• (V Xl w,X 2 ) + (X x , V X2 w) = 0, 

for all horizontal vector fields X\ and X 2 . 

Let v be a tangent vector of M. Then the projections of v onto 
the distribution D and its orthogonal complement D 1 - are called the 
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horizontal part arid the vertical part v VCT of v, respectively. Let 
/ : M — > K. be a smooth function. For the notation convenience, we 
also denote the horizontal part and vertical part of the gradient V/ by 
Vhor/ and V ver /, respectively. Let Vi,...,Vk be a frame in D which is 
orthonormal. Then the sub-Laplacian of / is defined by 

A hor / = ^(V. lIi V/,^). 

i 

Recall that the Ricci curvature Rc(u, v) is defined as the trace of 
the following operator w i— > (R,m(w, v)v, w). We define the horizontal 
Ricci curvature Rc hor by 

Rc ho >, v ) = ^ (Rm^, t>)^> 

i 

and the vertical Ricci curvature Rc vcr by 

Rc vcr (t>, v) = (Rm(«j, v)v, Ui) . 

i 

Let pt be a smooth positive solution of the following equation 

pt = A hoT p t + (VUi, VhorPt) + U 2 p t + Kp t \ogp t 

where U\, U2 are smooth functions on M and K is a constant. 

Let f t be the one-parameter family of smooth functions defined by 

f t = -2\og Pt -U l . 

A computation shows that f t satisfies the following equation 

(2.1) ft + ^|V hor / t | 2 = A hor / t + V + Kf t , 

where V = A hov Ui + KXJ X + ||V hor f/i| 2 - 2U 2 . 
We call a solution r of the problem 

r(t) = F(r(t)), r{t) 00 as t -> + 

stable if there is a family of solutions r e of the following 

f e (t) = F(r(t)) + e, r e (t) -> 00 as t ->■ + 

such that r e converges pointwise to r t . 

Finally, recall that a distribution is involutive if the Lie bracket of 
any two sections in the distribution is again in the distribution. The 
following is the main result of this paper. 

Theorem 2.1. Assume that the orthogonal complement D 1 - of the dis- 
tribution D is involutive and is given by the span of n — k horizontal 
isometries. Assume also that the following conditions hold: 

(1) Rc hor (v, v) + 3Rc ver (v hor , v hor ) > #iK or | 2 + I<2\v ver \ 2 , 



PAUL W.Y. LEE 



12 

or ; 



(2) Rc ver (v hor ,v hor ) <K 3 \v h , 

(3) V < K A , 

(4) A hor V + K 5 \V ver V\ 2 - 2K.V < K 6 , 

(5) a 3 (t) + ^>0, 

(6) ^(t) + + ( fll (t) + 1) (2K, - a 3 (t) - g|) = 0, 

(7) «fl - K 2 - 4K 3 + a 2 (t) (^g - ^ + ^ - f + = 0, 

(8) d 3 (t) + 222^1 + (X _ a3 (t)) ( a3 (t) + «j _ 2^) = 0. 

/or some positive constants Ki, K 2 , K 3 , K±, K5, Kq. 
Let r(-) be a stable solution of 

mtt t/ie condition r(t) — )■ 00 as £ — >■ + . 
Tnen 

A ft0 ,/ t (x) + ai (t)f t (x) + ^|V„ er / t (x)| 2 + a 3 (t)/ t (x) < r(t) 
for all t > and a// cc in M. 

In the case p t = Ah or pt, the proof of Theorem 12. II gives the following 
result which is one of the results in [3J. Note that, unlike [3J, the cur- 
vature conditions of the following result is written using a completely 
Riemannian notations. 

Corollary 2.2. Assume that the orthogonal complement D 1 - of the 
distribution D is involutive and is given by the span ofn — k horizontal 
isometries. Assume also that the following conditions hold: 

nfrorf„, „,-\ 1 nT>„ver{„, y ^ \ ^ |„, |2 

12 



(1) Rc hor (v, v) + 3i2c we >, or , v hor ) > K 2 \v verl 



(2) Rc ver (v hor ,v hor ) < K 3 \v h0 

for some positive constants K 2 , K 3 . Let p t be a smooth positive solution 
of the equation 

Pt = AharPt 

and let f t = —2\ogp t . Then 

A hor f t (x) + cf t (x)+2t (k 2 - |V, er / t (x)| 2 

4(c + ifkKl 



(cK 2 - AK 3 )(8K 3 - cK 2 )t 
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for all c > -^p, all t > 0, and all x in M. 

If we set c = in Corollary \2.2\ then the result further simplified 
to the following. 

Corollary 2.3. Suppose that the assumptions in Corollary \2.S\ hold. 
Then 

2 



A hor f t (x) + ^(*)+^|V_/ t (*)| 2 <^ 1 + ^ 



K 2 ' 3 1 - t 

for all t > and all x in M . 

In the case p t = A^ OT p t + (VUi, VhorPt) + U 2 p t , Theorem 12.11 gives 

Corollary 2.4. Assume that the orthogonal complement D L of the 
distribution D is involutive and is given by the span ofn — k horizontal 
isometries. Assume also that the following conditions hold: 

(1) Rc hor \v,v) + 3Rc ver (v hor ,v hor ) >K 2 \v ver \ 2 , 

(2) Rc ver (v hor ,v hor ) < K 3 \v hor \ 2 , 

(3) V < K A , 

(4) A hor V + K 5 \V ver V\ 2 <K 6 , 

y°) ZV5 cK 6 (8K 3 -cK 2 )> 

for some positive constants K 2 , K 3 , K A) K e and c > Let p t be a 

smooth positive solution of the equation 

p t = A hor p t + (VC/i, V 'horPt) + U 2 p t , 



ft = -2\og (H-Ui, c x = 2y/K 5 (K 2 - 4£a), andc 2 = ^^p. Then 

A hor f t (x) + cf t (x) + Cltan o h(C2t) |V^(x)| 2 < — coth(c 2 t) + cK 4 

2 C 2 

for all t > and all x in M. 

Next, we recall the definition of Sasakian manifolds and show that 
Theorem 11.31 is a consequence of Corollary 12.41 For a more detail dis- 
cussion of Sasakian manifolds, see [4]. 

Let M be a 2n + l dimensional manifold. A 1-form a on M is contact 
if da x is a non-degenerate 2-form on the kernel D of a (i.e. 

D x = {ve T x M\a{v) = 0} 

for each x). 

Let J be a (1, l)-tensor, w be a vector field, and a be a contact 
1-form on M. The triple (J, w, a) is an almost contact structure of M 
if the following conditions hold 

(2.2) J 2 (v) = -v, a(w) = 1, J(w) = 0, 
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where v is any vector in D. 

An almost contact structure (J, w, a) is normal if 

(2.3) [J, J](wi, w 2 ) + da(wi, w 2 )w = 

for any vector fields Wi and w 2 on M, where [J, J] denotes the Nijenhuis 
tensor defined by 

(2.4) [J, J](wi,w 2 ) = J 2 [wi, w 2 ] + [Jwi, Jw 2 }-J[Jwi, wi]—J[wi, Jw 2 }. 

An almost contact structure (J, w, a) together with a Riemannian 
metric (•, •) is called a almost contact metric structure if 

(Jvi, Jv 2 ) = (v 1 , v 2 ) and (w,w) = l, 

where V\ and v 2 are vectors in D. 

An almost contact metric structure is a contact metric structure 
if (wi,Jw 2 ) = da(uii,w 2 ). A Sasakian manifold is a manifold M 
equipped with a contact metric structure which is normal. 

An example of Sasakian manifolds is given by the Heisenberg group. 
The underlying space M of the Heisenberg group is the 2n + 1 dimen- 
sional Euclidean space M 2n+1 . In this case, the contact form a is given 
by 

1 n 

a = dz - - (yidxi - Xidyi) 

i=l 

where {x\, x n , yi, y n , z} are coordinates on M 2n+1 . 

In the Heisenberg group, the vector field w is given by w = d z . Let 
Xi = d Xi + \yid z and Yi = d Vi — \xid z . The Riemannian metric is 
defined such that {X±, ...,X n , Yi, ...,Y n , w} is an orthonormal frame. 
The tensor J is defined by 

JXi = -Yi, JYi = X i} and Jw = 0. 

Back to the general case, the Riemann curvature tensor Rm of a 
Sasakian manifold satisfies the following properties. 

Proposition 2.5. Let (M, J, w, a, (•, •)) be a Sasakian manifold. Then 
the fallowings hold: 

(1) (V wl J)w 2 = \ ({wi,w 2 ) w- (w,w 2 )wi), 

(2) V Wl w= -\Jw x , 

(3) w is a Killing vector field, 

(4) Rm(wi,w 2 )w= \({w,w 2 )wi - (w,Wi)w 2 ), 

for all tangent vectors Wi and w 2 on M. 
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A proof of the above proposition can be found in pj]. Note that the 
definition of the exterior differential da of a differential 2-form a used 
in this paper is 

da(w 1 ,w 2 ) = wi(a(w 2 )) - w 2 (a(w 1 )) - a([w 1 ,w 2 ])- 

This is different from the one in [4] and so the above formulas are also 
different from those in [1] by a multiplicative constant. 

Recall that the Tanaka connection V of a given almost contact metric 
manifold is given by 

1 1 1 

V Wl w 2 = V Wl w 2 + - (w, wi) Jw 2 + - (w, w 2 ) Jwi - - ( Ju>i, w 2 ) w. 



The corresponding curvature Rm is given by 



Rm(w u w 2 )w 3 = V Wl V W2 w 3 - V W2 V wl w 3 - V [wuW2] w 3 
and we denote by Rc the corresponding Ricci curvature 



Rc(^i,u;i) = trace(«>2 (llm(w 2 ,wi)wi,w 2 )). 

The following proposition shows that Sasakian manifolds provide 
examples to the main results. 

Proposition 2.6. The fallowings hold on a Sasakian manifold: 

(1) Rc ver (v,v) = Rc ver (v hor ,v hor ) = \\v hor \\ 

(2) Rc hor (v,v) + 3Rc ver (v hor ,v hor ) = Rc hor (v,v) + \\v hor \ 2 

"2 I ^wer I ~t~ R^iyhori^hor)' 

Proof. Clearly, we have 

Rc vcr (t>,i>) = (Rm(w, v)v, w) 

= (Rm(w,w hor )v hor , w) = Rc ver (> hor ,t; hor ). 
By Proposition 12.51 we also have 

(Rm(w,W h or)^hor, W) = (Rm(i7 hor , W)\V, V hor ) = ^hor| 2 - 

By Proposition 12.51 again, we have 

(Rm(wi 1 t))u ) w 2 ) = (Rm(wi,tw)tW,W2) 

+ (w, v) (Rm(ti)i, v)w, w 2 ) + (w, v) (Rm(w 1 , w)v hor , w 2 ) 

(2.5) _ (R m ( Wl]V ) V] W2 ) + _ ( Wj v f W2 ) 

K'vcrl 2 

= (Rm(Kii,i) lor )u hor , w 2 ) + V J (wx, w 2 ) . 
By the definition of Tanaka connection, 
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is horizontal for any vector fields X and Y. Therefore, by Proposition 

(V^V^^hor, Vj) = (V Vl V VhoT V hol ,Vj) 
= (V^V^hor - (V^Whor, W) w),Vj) 

(2.6) i 

= (V«, Vv hol Vhor, Vj) + - (V Vho V hor ,w) (JVi,Vj) 

Here we extend v to a vector field and still call it v. 
Similarly, we also have 

= (V^V^Vhor,^-) + - (V Vl V hor ,w) (Jv,Vj) 

By the definition of Tanaka connection, it also follows that 

(V[ Bi ,^]i;hor^i) 

(2.8) i 

= (VK.twFhor,^') + - (W, [Vi,V hoI ]) (Jv,Vj) . 

Therefore, by combining (12.6jl . ( I2.7jl . and f 12 .81) . we obtain 
(Rm^tw^hor,^-) = (Rm(vi,v hor )v hoi ,Vj) 

(2.9) + \( Jv > v i) ( Jv > v j) ~ 2 ( w ' t *' Vhor ]) u ^ 

3 

= (Rm(l) i , -y h or)fhor, Wj) + ^ Uj) ( Jt>, Uj) . 

Since VxY is horizontal, we also have 

(2.10) (Rm(w, v hor )v hor , w) = <Rm(w, v hoi )v hoT , w) = 

for any vector field w. 

Therefore, by (12. 5p . (12. 9p . and ( I2.10p . the second assertion follows. 

□ 

Proof of Theorem 11.31 It follows immediately from Corollary 12.31 and 
Proposition EE] with c = K 2 = § , K 3 = ±, and k = 2n. □ 

Proof of Corollary 1.4 Let 7(-) be a minimizer of the functional 
7(-)^mf£ 1 ^|7( S )| 2 + ^( 7 ( S ))^, 
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where / ranges over all smooth curves 7 (-) such that j(s) G for 
all t in [sq, Si]. 
By Theorem 11.31 

(i + j t M j(t)) = { 1 + l) /t(7(t)) + i 1 + 1) < Vto *fr(t)),7(t)> 

<-l\V ho J t \l + V(x) + ^coth(c 2 t) + ^+(l + -) (V hor /*(7(t)),7(t)). 
2 C2 n \ n J 

By Young's inequality, we have 

( 1 + -) ^/t(7W) < V{<y(t)) + *± coth(c 2 t) + ™± + 1 ( 1 + | 7 (t)| 2 . 
\ n J at C2 n 2 \ n J 

By integrating the above inequality, we obtain 

fl + -) (fsM) - fM) < ^ln (^^r) + (l + ^'^(xo.xO- 
V ny C2 V smh ( c 2 s o)/ V n / 

This result follows from this. □ 



3. Parallel adapted frames 

In this section, we define convenient adapted frames along a path 
called parallel adapted frames and use it to see how the linearization 
of a flow changes. 

Let M be a compact Riemannian manifold of dimension n equipped 
with a distribution D (i.e. a sub-bundle of the tangent bundle) of rank 
k. An orthonormal frame V\, Vk, «i, u n -k in the tangent space 
T X M at a point x is an adapted frame if v±, ...,Vk is contained in the 
space D x . 

Lemma 3.1. Let 7 : [0,T] — >■ M be a smooth path in M. Then there 
exists a 1-parameter family of adapted frames 

{vx(t), V k (t),Ui(t), u n _ k (t)} 

along 7 (t) such that 

(1) Vi(t) is in the orthogonal complement of D^u) f or 1 < * < k, 

(2) Uj(t) is in D^u) f or 1 < j < n — k. 

Here v denotes the covariant derivative ofv(-) along 7 (-). 

Moreover, if {ui (t) , . . . , u n -k (i) , v± (t) , . . . , v k (t) } is another such frame, 
then there are orthogonal matrices and of sizes k x k and 
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(n — k) x (n — k) ( independent of time t ), respectively, such that 

k 

v l (t) = Y,0 ( ll ) v l (t), l<i<k, 



1=1 

N-k 



s=l 

Let ux, ...,u n -k,vi, ...,Vk be an adapted frame in T 7 ( )M. According 
to the above lemma, there is a unique 1-parameter family of adapted 
frames U\ (t) , . . . , u n _k (t) , V\ (t) , . . . , v k (t) defined along 7 such that Vi (0) = 
V\ and Uj(0) = Uj for each % — 1, k and each j = 1, n — k. We call 
Ui(t), ...,u n -k(t),Vi(i), ...,Vk(t) a parallel adapted frame defined along 
7 which starts from Ui, u n _fe, i>i, Vf.- 

Remark 3.2. The notion of parallel adapted frame is a generalization 
of parallel transported frame in Riemannian geometry to the present 
setting (see the end of this section for detail). The construction of the 
above frame can be generalized to a more complicated setting where 
we have a smoothly varying flags of subspaces instead of just one dis- 
tribution. 

Let (p t be the flow of a vector field X t defined by ip t = X t (<pt) and 
<Pa(x) = x. The linearization dipt of (pt satisfies 

^j-d<p t (w) = S7dip t {w)X t . 

Here, we use 4 to denote covariant derivative along 1 1— >■ (pt(x). 

Therefore, the change in dip t is completely determined by the (1, 1)- 
tensor w H- V w X t . We will investigate the equation satisfied by this 
tensor. Let N(t) be the matrix representation of the (1, l)-tensor w 1— > 
V w X t with respect to the above parallel adapted frame at time t. More 
precisely, the r/'-th entry Nij(t) of N(t) is defined by 

[V Ui (t)X t , Uj(t)) if % < n - k and j <n - k 

[V Ul ( t )X t , Vj- n+k (t)) if % < n - k and j > n-k 

(V Wi _ n+fc (t)X t) Uj(t)) if z > n-k and j <n-k 

y v t _ n+k {t) x t-> Vj- n+k (t)) if % > n - k and j > n - k. 

Similarly, let R(t) and M(t) be the matrix representations of the 
bilinear form 

w (Rm(w, X t )X t ,w) 

and the (1, l)-tensor 

w^V w {X t + V Xt X t ), 



Nij(t) 
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respectively, with respect to the given parallel adapted frame at time 
t. Finally, let 

0( n -fc)x(n-fc) W(t) 

-W(t) T kxk 
where W(t) is the (n — k) x k matrix with r/'-th entry equal to 



Lemma 3.3. The 1 -parameter family of matrices N(t) satisfies the 
following matrix Riccati equation 

N(t) = -N(t) 2 - N(t)W(t) - W(t) T N(t) - R(t) + M{t). 

Finally, we split each of N(t), R(t), and M(t) into four pieces 

A/" 00 (t) Af i(t) \ r?u\ - ( n oo(t) n i(t) 



N{t) - I Mo(*) Mi(*) J ' m ~ \ K 10 (t) Knit) 



M{t) 



Moo(t) M 01 (t) 
M w (t) M n (t) 

where Aoo(i), ^oo(^), and A^oo(^) are of size (n — k) x {n — k). 

The following, which will be used in the later sections, is an imme- 
diate consequence of Lemma 13.31 

Lemma 3.4. The 1-parameter family of matrices N(t) satisfies the 
following 

Jif n (t) = -M u (t) 2 - Af 10 (t)Af i(t) 

-Af 10 (t)W(t) - W(t) T Af 01 (t) - n u (t) + M n (t). 

Proof of Lemma [3J[ Let (vi(t), Vk(t),Ui(t), u n -k(t)) be a 1-parameter 
family of adapted frames along the path j(t). Let 

(vx(t), ...,W fc (t),tii(t), ...,u n - k (t)) 

be any other such family. Then 

k 

Vi(t) = ^20 a (t)vi(t), l<i<k, 
i=i 

where 0(t) are orthogonal matrices of size k x k. 

Let A(t) be the k x k matrix with ij-th entry equal to (vi(t), Vj(t)). 
Then (vi(t), vj(t)) = for each 1 < i,j < k if and only if 

0{t) = -0(t)A{t). 

Since A(t) is skew symmetric, we have a solution 0(t) to the above 
ODE. Moreover, any solution is determined by its initial condition. 
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This proves the result for «»(•). A similar procedure gives the result for 

Ui(-). □ 

Proof of Lemma \3.3l Let 

V(t) = {u 1 (t),...,u n - k {t),v 1 (t),...,v k {t)) T 

and 

V(t) = (Ui(t), ...,Un- k (t),Vi(t), ...,i) k (t)) T , 

where {ui(t), ...,u n - k (t),Vi(i), ...,Vk(t)} is a parallel adapted frame de- 
fined along the curve t H- ipt( x ) as i n Lemma [37TT 
By Lemma [3. 1\ we have 

If we differentiate the above equation once more, then we obtain 

v(t) = ms(t)v(t) + w(t)v(t) = (ms(t) + w(t) 2 )v(t). 

Let = (d(ft(v o (0)) , dip t (v2n{0))) T and let A(£) be the matrices 
defined by 

$ t = A(t)V(t). 

It follows that 

(3.1) j$ t = (A(t)+A(t)mw(t) 

and 

(3.2) = + 2A(t)W(t) + A(t)SIJ(t) + A(t)2H(t) 2 )y(t). 

On the other hand, if we let 7(5) be a path such that 7'(0) = 1^(0), 
then 

D D 2n 

(3.3) -d<p t ( Vi (0)) = -X t (^( 7 ( S ))) _ n = X)^(t)V sjW X t (^(x)). 



s=0 

j=0 



By the definition of Rm, we also have 

--^ t (^(0)) + ^A J (t)Rm(^(t),^(^(x)))X t (^(x)) 

j=0 

(3-4) =^§^(^(7(,)))| 

as at \s=o 

In 

= J2 A dt) (v,,(*)(X t + V* t X t )) fo(x)). 

j=0 
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By (13. ip and ( 13.31) . we have 

(3.5) N(t) = A(t)- l A(t) + 2H(t). 

By (13. 2 p and (13. 4p . we also have 

A{ty 1 A{t) + 2A{t)- 1 A{t)W{t) 

+ 23J(t) + 2U(t) 2 + R(t) - M{t) = 0. 

If we combine ( 13 .5D and ( 13. 6p . then we obtain 

JV(t) = -N(t) 2 - N(t)W(t) - W(t) T N(t) - R(t) + M(t) 

as claimed. □ 

Before ending this section, let us discuss the relationships between 
parallel transported frames and the Tanaka connection. In the usual 
Riemannian case, if v\(t), ...,v n (t) is a parallel orthonormal frame de- 
fined along a path 7, then one can define the covariant derivative v(t) 
of a vector field 

v(t) = ai(t)vi(t) + ... + a n (t)v n (t) 

defined along j(t) by 

v(t) = ai(t)ui(t) + ... + d n (t)v n (t). 

It is, of course, well-known that the covariant derivative is closely re- 
lated to the corresponding Levi-Civita connection. 

Similarly, one can define certain covariant derivative corresponding 
to the above parallel transported frames. More precisely, let 

w(7(*)),«i(*)> ->f2»(*) 

be a parallel transported frame defined along a path 7 in an almost 
contact metric manifold. If 



v{t) = ao(*)w(7(*)) + ai(t)vi(t) + ... + a 2n (t)w 2 n(£) 



is a vector field defined along 7(£), then the covariant derivative -j. cor- 
responding to the parallel transported frames of v(t) along 7 is defined 
by 

— V (t) = Oo(t)w(7(t)) + ai(*)Ul(t) + ••• + 02n(*)«2n(*)- 

ctr 

Note that the definition of £ is well-defined. 

The following lemma gives some basic properties of % and some of 
its relationships with the Tanaka connection V. Since it is not needed 
for the rest of the paper, the proof is omitted. 
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Lemma 3.5. Let w(t), Wi(t), and w 2 (t) be three vector fields defined 
along a curve 7 in an almost contact metric manifold. Let c\ and C2 be 
two constants and let a(t) be a smooth function. Then 

(1) £ (ciwi(t) + c 2 w 2 (t)) = ci|wi(t) + c 2 gw 2 (t), 

(2) f t a(t)w(t) = a(t)w(t) + a(t)§w(t), 

(3) i«=0, 

(4) i/ &o£/i Wi(t) and w 2 (t) are contained in D, then 
^( Wl (t),w 2 (t)) = (^ Wl (t),w 2 (t)\ + /wi(t),^w 2 (t) 

(5) if Y is a vector field contained in D, then 

^Y^it)) = V l{t) Y - l -{w^(t)) JY- 



4. Distributions with Transversal Symmetries 

In this section, we assume that the orthogonal complement of the 
given distribution D is spanned by n — k horizontal isometries denoted 
by Wi,...,w n _fe We will also assume that the vector field X t in the 
previous section is the horizontal gradient Vhor/t of a one-parameter 
family of functions f t defined on the manifold M and specialize Lemma 
13.41 to this case. 

Lemma 4.1. Assume that the orthogonal complement D 1 - of the dis- 
tribution D is involutive and is given by the span of n — k horizontal 
isometries. If ift is the flow of a time- dependent vector field Vhorft- 
Then 

d X f 1 

—A hor f t (ipt) < —^(A hor f t (^ t {x))) 2 + A hor (f t + -\V hor f\ 2 

- Rc h ; t { x) (v f t , vf t ) + R C ; e t r {x) (v hor f t , v hor f t ) 

k 

Proof. Let us use the notation as in Lemma 13.31 with X t = Vhorft- 
Recall that ui(t), ...,u n -k(t),vi(t), ...,Vk(t) is a parallel adapted frame 
along the path t 1— > <^>t{ x )i where (ft is the flow of Vhor/t- By (1) of 
Lemma [6.11 we have 

tr(MiW) = J2( V Mt)^orft,Vi{t)) 

(4.1) 
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By (1) of Lemma [6J] and the symmetry of the Hessian, we also have 

tr(AT n (t) 2 ) = J2( V n(t)Vhorf,Vj(t)) (V Vj(t) V ho J, Vi (t)) 
hi 

= E « V ^) V M(*)) - <v WiW v ver /,^(t)» 

hi 

■ «V„ jW V/,^(t)> - (V„ iW V ver /,^(t)» 
= E (< V ^) V MW) + (V„. (t )V ver /,^(t)» 

• ((V w(t )V/, - (V % . (t) V ver /,^))) 

= E < v ^wv/,^(t)> 2 - E <v, jW v vcr /,^(t)> 2 . 

Therefore, by (11) of Lemma [6.11 and the Cauchy-Schwartz inequal- 
ity, we have 

(4.2) tr(A/- n (t) 2 ) > i(A hor / t (^(x)) 2 - Rc£fo(V w /, V vcr /). 
Let O(i) be a family of orthogonal matrices such that 

j 

It follows from (1) of Lemma [6. II that 
M^-(t) = (iii(t),Vj(t)) = E°*(*) (Vv hor/i w fc (^(x)),^(t)) 

= -E° ifc ^ (V^(t)W fc (^ t (x)), V hot Mvt{x))) 

k 

= E 0ifc W (wfeC^*^))' V Wj .(t)Vhor/t) = (V„ i(t )Vhor/t 5 Wi(t)) . 

Therefore, we have 

(4.3) W(t) T = Mo(*)- 

Therefore, the ij'-th component of iy(t) T A/oi(t) = A/io(£)A/oi(£) is 
given by 

E(^(*) ) ^(*))(V U! (t)Vhor/t ) ^(t)) 
I 

h rft W k{<Pt(x)),Vi(t)) Ol s (t) (Vw s Vhor/i(¥>t(a;)),Uj(t)) 

k,l,s 

= E ( v v h «/t w fc(¥'t( a; ))> «<(*)> (Vw fe Vhor/t(^(x)),^(i)) 
k 
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and so 



tr(Mo(tyVoi(t)) = tr(^(t) J JV i(*)) 

f\7 V7. AY \ 



( ^ 4 ' 4 ' ) = ^ ((V Vhor / t W fc ) hor , (V Wfc V h or/t)hor) 



A: 

The zj-th component of Afio(i)Nio(t) T is 

V v hoI /w fe > O ls (t) (V Vhor/t w s ,^(t)) 

fc . s 

= Yl V Vhor/ W fe> (Vv hor / t Wfe, ^(t)) . 

Therefore, 

|Mo(*)| 2 = tr(MoWMo(*) T ) 

= K V V hol ./iWfc)hor| 2 t ( x ) = Rc™ r (x)(Vhor/t, Vhor/t)- 

k 

By (1) and (3) of Lemma [6.11 it follows that 



trCMn(t)) = (v Ui(i) (V hor / t + V Vhor/t V har /t),«iW < 

= A hor U + i|V hor/ | 2 ) - 2 <V Ui(t) V Vhor /V ver /, . 

On the other hand, we have, by (8) and (12) of Lemma [6.11 
J] <V„ i(t) V Vhor / t V ver / t , Ui(t)> = ( V n(t) ( w «/ v v hor / t w z ) , Vi(t)) 

i i,l 

= R<(.)(V/ t , v vcr /,) + £ (V hor (wJ t ), Vv te/( w,> wW 

= RcJ^V/t, Vvcr/t) + ((Vw ; V hor / t - V Vhor/t wOho r , Vv^/.W;)^^ . 

I 

Therefore, by combining this with OH]), flO]) . ( [Qjl . and (113]) . 

we obtain 

^A hor / t (^) = ^tr(Mi(t)) 

= -tr(AT n (t) 2 ) -2tr(Mo(t)A^ i(t)) - |Mo(*)| 2 - tr(TZ n (t)) + tr(M n (t)) 
< -I(A hor / 4 (^(x))) 2 - Rc^ s) (V/ t) V/ t ) + RC (x) (V hor / 4 , V hor / t ) 

V hor / t Wfe)hor, (V Wfc V] a or/t)hor) v , t ( a; ) + A hor ^/f + -|V h or/| 2 ^ • 
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□ 



5. Proof of the main results 

This section is devoted to the proofs of the main results. We begin 
with 

Lemma 5.1. Assume that the orthogonal complement D 1 - of the dis- 
tribution D is involutive and is given by the span of n — k horizontal 
isometries. Then the followings hold: 

(1) |(/ t (p t )) = ~fM) + 2A hor f t (^ t ) + 2V(<p t ) + 2Kf t (<pt), 

(2) i(Mn)) = ^ftM + KftM, 

(3) i (||V„ er /4J = A hor (\\V ver f t \ 2 ) (<p t ) 

~ Et \(V Wi Vhorft - ^V hor f t Wi) hor \l t + (V ver f t ,VV) Vt . 

Proof of Lemma \5. 1[ The first assertion follows from (12. ip . By (I2.ip . 
we have 



d 
dt 



iftM) = j t f-^|V hor /4 t + A hor / t (^) + Vfa) + KM^yj 

- ^Vhor/t, Vhor/t^) - - (V| Vhor/t| 2 , Vhor/t)^ + A hor f t {ft) 



ft 2 

2 

tpi 



+ (V hor A hor / t , V hor / t ) + (VV, V hor /) + Kf t (<p t ) + K\ Vhor/, 
= A ho J t (ip t ) + Kf t (<p t ). 

By (1) and (8) of Lemma [6.11 we have 

(V hor (wi/), V hor /) = (V Wi V hor / - Vy^/Wj, V hor /) 

= ^vQ|V hor /| 2 ),w^. 
Therefore, by combining this with (12. ip . we have 
j t Q|Vver/t|Jj =^2wJ t (ip t )j t (w i f t (ip t )) 

= ^WiMVt) CwiftM + (Vhor/t, Vhor(Wi/t)) 



^Wi/t^w^A/* + V + Kf t )( Vt ) 

i 

(V vei f t , VA hor f t + VV) + K\V vcr f\l t . 
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Finally, by (10) of Lemma [6. 1[ we obtain 



\V vei f t \l t = A hor -|V vcr / 4 | 2 {ipt) + K\V vc J t \l t 



□ 



Proof of Theorem \2.1[ By Lemma [4.11 and (12) of Lemma [6.11 



d A«/tfo) < " ^A hor / t (^) + A hor (f t + hv hOT ff 



at k k \ 2 

- Rc^CV/t, V/ t ) - 3Rc; c /(V hor / t , V hor / t ) 



4 5Z ( v v hor / t w fc , V Wfc V hor / t - V Vhor / t w fe )^ . 



By Young's inequality and (12) of Lemma [6.1[ the above inequality 
becomes 



|A hor /^ t ) < ^f- - ^A hor /^) + A hor U + i| V hor jf) M 
+ a 2 (t) V |(V Vhor/t w fc - V Wfc V hor / t )hor|J t + -^-Rc™{V hoi f t , V hor f t ) 
- RcJ° r (V f t , Vf t ) - 3Rc; c t r (V ho J 4 , V h0I f t ). 



Using the first and the second assumptions, 



d a r/..N^04(*) 2 2a 4 (t) A j? /,„ \ i a A 1 1T7 



— A ho Jt(tp t ) < — —A ho Jt(tp t ) + A hor ( / t + -| Vhor/| ) (Vt) 



+ a 2 {t) ^ l(Vv ho ,f t Wfc - V Wfc Vhor/t)hor 



+ ^|v hor /4 t - #i| WtlJ, - ^ 2 |v vcr / t |J t . 
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Let F t (x) = A hoT f t (x) + ai (t)f t (x) + ^f\V vc J t (x)\ 2 + a 3 (t)f t (x). By 
Lemma [5. 11 we have 

j/ t fa) < ^f- ~ 2 -^A ho J t fa) + A hor U + ^|V hor /| 2 ) fa) 
+ ^r\V ho J t \l t - K^Ml - K 2 \V vc J t \l + a^Ufa) + ^-\V vc J t \ 2 



+ a 3 {t)f t fa) + ai(*)A hor / t (^) + a x (t)Kf t fa) + a 2 (t)A hor Q|V ver ./t| 2 ^ (p t ) 

+ a 2 (t) (Vver/t, VV)^ + a 3 {t){-f t fa) + 2A bat f t (<p t ) + 2F(^) + 2Kf t (<p t )). 
By (12. ip and collecting terms, we obtain 

j/tfa) < + A hor F t + A hor V + 2a 3 (t)Vfa) 



2 



+ ( a3(t) _^ +A - +2 ( S| _, 1 )) Wl 

+ ( di(t) + ai(*)# - a 3 (t) - 2 ( - ) f t fa) 



[ a 3 {t) + 2Aa 3 (t) + 2K ( - K x \ \ f t fa) 



a 2 (t) 
i_ 

\a 2 (t) 
+ a 2 {t) (V vc Jt,VV) Vt . 

By Young's inequality and the definition of F t , we obtain 

^-F t fa) < + A hor F t (^) + A hoi Vfa) + 2a 3 (t)Vfa) 

at k 

+ K 5 \ V vcr V| I + (a 3 {t) - 2 -^ + K+^- 2K^ F t fa) 
Vfa) + (^-K 2 - 4K 3 ) |V ver / t | 2 



J- a (A ( ^ -u ^ K +lc\\K7 f |2 

+ a2(t) \m --t + — -y + Ki ) |Vvcr/ ^ 

. i{t) + 2^^) + + 1} ^ _ ^ _ ^ 

+ f d 3 (t) + + (A - 03(f)) (a 3 (t) + ™jL - 2A,) ) /^). 
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By assumptions (6), (7), and (8), the inequality becomes 



2 



2 

ft 



,M<p t ) < -^y- + A h0I F t (tp t ) + A hor V(tp t ) + K 5 \V vcr V 
at k 

( , \ 2aJt) 8K 3 *\ , . 

+ 2 («•<«> + -*)"<<*>■ 

By assumptions (3), (4), and (5), 

TMft) < ^f- + ^oM^t) + K 6 + 2 (a 3 (t) + ^r) K, 
at k \ 0,2{t) J 

Let e > and let r e (-) be a solution of 



fe(t) = -^r~ + K 6 + 2 a 3 (t) H rf- 

k V a 2(t)/ 

a 3 (t)- 2 ^ + K+^-2K 1 ) r e (t) + e. 

with condition r e (£) — >• oo as £ — > 0. 

Let to > be the first time where there is a point x in M satisfying 
F t (¥t 0)) = r(t )- Then 

fe(t ) < ^Ftfa) < + K 6 + 2 (a 3 (t ) + K, 

at k \ a 2 {t ) J 

a 3 (t ) - ^ + K + _ 2Kl ) r e (t ) < r e (to) 

k a 2 (£ ) / 

which is a contradiction. 

Therefore, F t (x) < r e (t) for all t > and all x in M. The result 
follows from stability of r. □ 

Proof of Corollary \2.2\ This follows from Theorem 12. II by setting K\ 

c+l)kh 

{cK 2 -4:K3)(_8K3-cK 2 )f 

Proof of Corollary \2.3[ If we set K\ = K — 0, a 3 = 0, a x = c is a con- 
stant, a 2 (t) = citanh(c 2 t), a 4 (£) = **5g|p, and K 5 = ^skI-^L) 
in Theorem 12.11 Then a computation shows that r(t) = ^ coth(c 2 t) + 
CK4 and the result follows from Theorem 12.11 □ 



K = 0, a 3 = 0, ax = c, a 4 (t) = ±(f ±iM3, a 2 (t) = 2 (AT 2 - *fa) £, and 

' VV fcJCo— lA < U8A<-cAc> > >r 
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6. Appendix 

In this appendix, we provide the detail calculations that we used in 
the proof of Theorem 12 . 1 1 



Lemma 6.1. Suppose that the assumptions in Lemma 4jJ_ hold. Let 
W\, w n -k be the horizontal isometries which span D- 1 . Let X\ and 
X 2 be vector fields contained in D and let Z be a vector field contained 
in D 1 - . Then the followings hold: 

(V Xl Z,X 2 ) = -(V Xa Z,X 1 ), 
V Vi Vi is horizontal, 
V Wi Wj is vertical, 
W Vj Wi is horizontal, 

Vv hor fVhorf = |V 'fc or |V 'horf\ 2 ~ 2(Vy hor /V verf) hor, 

E j (V[v k ,v J ]Vf,v j ) = o, 

Ej (^[Wi^horf, Vj) = ~ J2j (V % VlJ, [W h Vj\), 
V h or{Wif) = (V^Vfeor/ - Vv horf Wi) hor , 
A hor (Wif) = Wi(A hor f), 

A hor {\\v ver f\ 2 ) = (v ver f,vA hor f}+^j(v w y hor f-v Vhorf w, 

Rc (V ' verj ') ^ verf) \\^Vj^verf)hor\ > 

Rc hor (Vf, «%) = E, (V Ui V VhorfWi , v 3 ), 
(Rmiv^w^Wj^j) = (V Vj Wi, V Vj Wj). 



(1 
(2 
(3 
(4 
(5 
(6 
(7 
(8 
(9 
(10 

(11 
(12 
(13 



i j hor | ; 



Proof. By assumption, we have 

(V Wl X 1 ,X 2 ) + (V Wl X 2 ,X 1 ) = w, (X 1 ,X 2 ) = ([w i ,X 1 ],X 2 )+([w i ,X 2 ],X 1 ) 

This gives (1). 

It follows from (1) that 

(V Vi Vi, Wj) = - (vi, V Vi Wj) = 

which is (2). 

Since D is involutive, we also have 

(Vw.Wj.Dfc) = - (Wj, V Wi V k ) = - (Wj, V Vk Wi) 

= (V Vk Wj,Wi) = (V Wj U fc , Wj) = - (l' k , V Wj Wi) = - (Dfc.Vw.Wj) . 

This gives (3). 
Since 

(V^.Wi, w fc ) = (V Wi Vj,Wk) = - (f„V Wt w/i) = 0, 

(4) holds. 
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The statement (5) follows from 
UV|V hor /| 2 ,Xx\ = (V Xl V hor f,V ho J) 

= (v Xl v/, V hor /) - (v Xl v ver /,v hor /) 

= (V Vhor /V hor /, X x ) + 2 (V Vhor /V ver J, Xt) . 
Since Wj is a horizontal isometry, [wj,^-] is horizontal. Therefore, 
E (V[w ijWj -]V/, = J] (V«, V/, [w <9 Wi ]> 

= E V/, Ufc) (ffc, [wi, u 3 -]> 

= - E ( v «j V /, «fc) [w 4 , ^]) = 

which is (6). 

It follows from (1) and (6) that 

= E < v [wi«]V/, «i) = E < v ^ t w - "iO 

= E ( V ^'j V hor/, [Wi, U,-]) + (V^ V vcr /, [W 4 , U,-]) 
j 

By (6), we also have 

= E (^[wi^jVhor/, ^) + E ( V 
j j 

Therefore, (7) follows. 

Statement (8) follows from 

(V( Wi /),^) = (V Wi V/,Wi) + (V hor /, V Wi Wi) + (V ver /, V Vj Wi) 

= (V Wi Vf,Vj) + <V hor /, V Vj Wi) + (V^ 
= (V Wi Vhor/,Vj) - (Vv hM /Wi,%> • 
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By (1) and (8), we have 

A h or(Wi/) = ^ <V^(V Wl V hor / - V Vhor /Wi), Vj) 

3 

= Yl ( V w, V,, 3 Vhor/, Vj) - (Vy^/V^.Wi, Vj) 
3 

+ Yl ( V K,w i ]V h or/ - V h , Vior /]Wi, U,-) . 
3 

By (1), (3), and (4), the above equation becomes 

A hor (Wi/) = ^(V Wi V^.Vhor/,^) - (Vvhor/V^.W;,^) 
3 

+ ( v h,w I ] V hor /, vj) + (V^w i; V^Vhor/) - (V^.Wi, V Vhor /%) • 
By (1) and (7), 

j 

3 

= < v w, v^. V hor /, uj) + ^ (v„.v hor /, V Wi «j-) . 

j j 

This, together with (6), gives (9). 
By (8) and (9), we have 

A hor Q|V ver /| 2 ) = J2(viWjf) 2 • ^w,/-A,,„ivv / /) 

= |V hor ( Wj /)| 2 + (V ver /, VA hor /> 

3 

= K V w 4 Vhor/ - Vv^/wOhorl 2 + (V w /, VA hor /> 

which is (10). 

By (3), Vv V er/^ver/ is vertical. Therefore, by (1), we have 

Rc h ° r (V vcr /, Vver/) = - <Vv vcr /V. t , Vvcr/, ^) " ^ ( V ^W] Vyer/, • 

J 3 
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By (1), the above becomes 

Rc h ° r (V vcr /, V vcr /) = ((V % V ver /, V Vvcr /^> + (V^.Vver/, [Vj, V vcr /]hor» 

i 

= ^ (( V ^' V vcr/, (V Vvcr /^')ver) + | ( V Vj V vcr /)hor | 2 ) • 
3 

It follows from (4) that Vv vcr /t)j is horizontal. Therefore, (11) holds. 
By (1) and (3), 

Rc hor (V/, Wi ) = (V^Vv^/w,,^-) + (V^.Wi, V vf Vj) - <V hiV /]W 4 ,^> 

3 

= J2 < v ^ v Vho r / w ^ v j) + (V^.Wi, V V /^> + (V^.Wi, [«,-, V/]) 
j 

= E < V ^ V v hor /w i , ^> + <V, 3 w„ V VJ V/) . 

3 

Statement (12) follows since 

(V„. Wij V„ 3 . V/> = E (V^Wi, « fc ) (v fc , V, 3 V/) = 

by (1). 

By (1) and (3), we have 

(Rm(vj, Wi)w k ,Vj) = (V Wj .wjfe, V Wi v,-) + <V^.w fc , [t^wj) 

= (V^Wi, V Ui w fc ) 

which is (13). □ 

References 

[1] B. Andrews: Harnack inequalities for evolving hypcrsurfaccs. Math. Z. 217 

(1994), no. 2, 179-197. 
[2] D. Bakry, M. Ledoux: A logarithmic Sobolcv form of the Li-Yau parabolic 

inequality. Rev. Mat. Iberoam. 22 (2006), no. 2, 683-702. 
[3] F. Baudoin, N. Garofalo: Curvature-dimension inequalities and Ricci lower 

bounds for sub-Riemannian manifolds with transverse symmetries, arXiv: 

1101.3590 (2012). 

[4] D.E. Blair: Ricmannian geometry of contact and symplectic manifolds. Second 
edition. Progress in Mathematics, 203. Birkhuscr Boston, Inc., Boston, MA, 
2010. 

[5] H.-D. Cao: On Harnack's inequalities for the Kahlcr-Ricci flow. Invent. Math. 

109 (1992), no. 2, 247-263. 
[6] H.-D. Cao, S.-T. Yau: Gradient estimates, Harnack inequalities and estimates 

for heat kernels of the sum of squares of vector fields. Math. Z. 211 (1992), no. 

3, 485-504. 



DIFFERENTIAL HARNACK INEQUALITIES ON SASAKIAN MANIFOLDS 27 

[7] B. Chow: On Harnack's inequality and entropy for the Gaussian curvature 

flow. Comm. Pure Appl. Math. 44 (1991), no. 4, 469-483. 
[8] B. Chow: The Yamabe flow on locally conformally flat manifolds with positive 

Ricci curvature. Comm. Pure Appl. Math. 45 (1992), no. 8, 1003-1014. 
[9] R. Hamilton: A matrix Harnack estimate for the heat equation. Comm. Anal. 
Gcom. 1 (1993), no. 1, 113-126. 
[10] R. Hamilton: The Harnack estimate for the Ricci flow. J. Differential Geom. 

37 (1993), no. 1, 225-243. 
[11] R. Hamilton: Harnack estimate for the mean curvature flow. J. Differential 

Geom. 41 (1995), no. 1, 215-226. 
[12] P.W.Y. Lee: On the Li-Yau estimate, arXiv: 1211.5559 (2012). 
[13] P. Li, S.-T. Yau: On the parabolic kernel of the Schrdinger operator. Acta 

Math. 156 (1986), no. 3-4, 153-201. 
[14] J. Moscr: A Harnack Inequality for Parabolic Differential Equations. Commun. 

Pure Appl. Math 17 (1964), 101-134. 
[15] L. Ni: Monotonicity and Li-Yau-Hamilton inequalities. Surveys in differential 
geometry. Vol. XII. Geometric flows, 251-301, Surv. Differ. Geom., 12, Int. 
Press, Somcrville, MA, 2008. 
E-mail address: wylee@math.cuhk.edu.hk 

Room 216, Lady Shaw Building, The Chinese University of Hong 
Kong, Shatin, Hong Kong 



